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Abstract: The scale of neutrino masses and the Planck scale are separated by more 
than twenty-seven order of magnitudes. However, they can be linked by imposing 
the stability of the electroweak (EW) vacuum. The crucial ingredient is provided by 
the generation of neutrino masses via a seesaw mechanism triggered by Yukawa inter¬ 
actions between the standard model (SM) Higgs and lepton doublets and additional 
heavy right-handed neutrinos. These neutrinos participate to the renormalization 
group (RG) running of the dimensionless SM couplings, affecting their high-energy 
behavior. The Higgs quartic coupling is dragged towards negative values, thus alter¬ 
ing the stability of the EW vacuum. In the usual type-I seesaw model, this effect is 
too small to be a threat since, in order to comply with low-energy neutrino data, one 
is forced to consider either too small Yukawa couplings or too heavy right-handed 
neutrinos. In this paper we explore this general idea in the context of low-scale 
seesaw models. These models are characterized by sizable Yukawa couplings and 
right-handed neutrinos with mass of the order of the EW scale, thus maximizing 
their impact on the RG flow. As a general result, we find that Yukawa couplings 
such that Tr (YjY„) > 0.4 are excluded. We discuss the impact of this bound on 
several observables, with a special focus on the lepton flavor violating process /i -»• ey 
and the neutrino-less double beta decay. 
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1 Introduction 

The discovery of the Higgs boson at the LHC [1, 2] represents a cornerstone in 
particle physics, finally achieved after a nearly half-century quest. Paradoxically 
though it may seem, this astonishing discovery did not clarify the true origin of the 
electroweak (EW) symmetry breaking: after LHC run I, any clue of new physics 
beyond the Standard Model (SM) is still missing. 

To crown it all, the measured value of the Higgs boson mass, M h = 125.09 ± 
0.21 stat . ± 0. ll S yst. GeV according to the latest combination of both ATLAS and CMS 
experiments [3], is consistent, in spite of any naturalness argument, even if the SM is 
extrapolated up to the Planck scale. This remarkable conclusion is based on the fact 
that the value of the Higgs boson mass is caught in the vice of two solid theoretical 
constraints if one tries to extend the SM up to arbitrarily high energies by means 
of its Renormalization Group Equations (RGEs). On the one hand, for large values 
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of M h , the Higgs quartic coupling A is driven towards greater and greater values, 
eventually becoming non-perturbative. On the other one, for small values of M h , 
A is driven towards smaller values, eventually becoming negative. The occurrence 
of the latter condition may have dramatic consequences since it implies that the 
EW vacuum is only a local minimum, and there exists a global minimum at higher 
energy in the Higgs potential [4-17]. For a sufficiently old Universe, nothing would 
prevent the vacuum expectation value (vev) of the Higgs field to be at the global 
minimum, as a mere consequence of a tunneling transition between the two vacua [18- 
20]. Fortunately, Nature spared us from this unspeakable catastrophe: in the SM 
with M h ^ 125 GeV A does become negative at high energies, but its absolute value 
remains far too small to trigger the tunneling. In other words, in the SM the EW 
vacuum is metastable [21-23]. 

The metastability of the EW vacuum holds under the assumption that the SM is 
the full theory all the way up to the Planck scale. This argument poses a legitimate 
question: what is the role of new physics beyond the SM with respect to the stability 
of the EW vacuum [24-35]? First and foremost, to address this question we have 
to carefully specify which kind of new physics we are referring to. Undeniably, the 
most compelling case is provided by the observation of non-zero neutrino masses 
and oscillations. This fact cries out for the inclusion in the SM of a mechanism 
responsible for the generation of neutrino masses [36]. A very economical way is 
provided by the so-called type-I seesaw mechanism [37-40]. In a nutshell, the SM is 
enlarged with the addition of three heavy right-handed neutrinos 1 with a Majorana 
mass Mr breaking lepton number, coupled to the SM Higgs doublet H and leptonic 
doublets L via a Yukawa coupling Y v . Integrating out the heavy fields, the effective 
light neutrino mass m u « v 2 Y u /Mr is generated after EW symmetry breaking, where 
v - 246 GeV is the EW vev. 

In principle, this additional Yukawa interaction pushes the SM towards the in¬ 
stability region, since, at least at one-loop, it exactly mimics the effect of the top 
Yukawa coupling on the running of A: it introduces a negative contribution which 
drags the Higgs quartic coupling towards negative values [22, 42-44], However, in 
order to have a sizable effect, the Yukawa coupling Y v must be of order one. This is 
a necessary but not a sufficient requirement. Indeed, it implies Mr ~ 0( 10 15 ) GeV 
if we assume m v ^ 0.1 eV for the light neutrino mass scale. This means that the 
right-handed neutrinos actively participate to the running of A only for values of the 
RG scale larger than Mr. Or, to put it another way, there is not enough time, in 
terms of RG evolution, to sizably alter the SM picture. This result does not change 
trying to lower the mass scale Mr since in this case, in order to reproduce the correct 
order of magnitude for the mass scale of light neutrinos, one is forced to consider 
Y v ~ (9(10- 5 ). Of course, this value is too small to modify the running of A. 

x An even simpler setup would imply only two sterile neutrinos [41], 
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However, this is anything but the end of the story. Despite its simplicity, the 
minimal type-I seesaw suffers from a penalizing phenomenological issue: it is not 
testable. As discussed before, it introduces either unaccessible large mass scale or 
extremely tiny Yukawa couplings. The necessity to overcome this unpleasant problem 
led to the introduction of extended seesaw models with both TeV-scale right-handed 
neutrinos and sizable Yukawa couplings [45-50]. The key feature of these models is 
that lepton number is softly broken via the introduction of extra singlet fermions in 
addition to the usual right-handed neutrinos. The mass Mr can be brought down to 
the EW scale without neither causing problem with low-energy neutrino phenomenol¬ 
ogy nor lowering the Yukawa coupling Y v . As a byproduct of this construction, a 
very rich low-energy phenomenology emerges. Potentially interesting signals include, 
for instance, lepton flavor violating radiative decays, deviations from EW precision 
observables, and production at colliders of heavy Majorana fermions. 

As already clear from these introductory comments, low-scale seesaw models 
may alter the metastability of the EW vacuum since they feature, at the same time, 
sizable Yukawa couplings and relatively low mass thresholds [51]. In this work we 
study in detail the stability of the EW vacuum considering the so-called inverse 
seesaw model (ISS) [45, 46] as a reference model with low-scale seesaw. We include 
in our analysis the constraints coming from low-energy neutrino phenomenology with 
the aim to provide a complete and realistic description of the physics involved. 

Schematically, this paper is organized as follows. In section 2 we discuss the 
most relevant features of the ISS model, with a special focus on the phenomenological 
constraints included in our analysis. In section 3 we study the stability of the EW 
vacuum in the ISS model, using the results of section 2. In section 4 we present our 
results, and finally we conclude in section 5. In appendix A we generalize our results 
to the linear and double seesaw models. 


2 The inverse seesaw model 


In this section we start our analysis describing the most relevant prerogatives of the 
ISS model. In section 2.1 we introduce the Lagrangian, the mass matrix and the 
corresponding diagonalization. In section 2.2 we discuss the constraints coming from 
low-energy neutrino phenomenology included in our analysis while in section 2.3 we 
present our strategy for the numerical analysis. 


2.1 The model 


In the ISS model the SM field content is extended to incorporate ur right-handed 
neutrinos N l R and riy singlet fermionic fields S' 1 . The Lagrangian in the ISS model is 
given by 


Ass - iNrY^Nr) + iSY(dY) 


NrY v H ] L + NrMrS + ^S c nsS + h.c. 


( 2 . 1 ) 
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where L = (L e ,L^,L r ) T represents the left-handed lepton doublets with the usual 
contents L l=e ' /1 ’ T - (tsil,Il) t while H is the Higgs field, with H = ■ Y v is the 

hr x 3 Yukawa matrix mediating the interactions between the SM leptons and the 
right-handed neutrinos while Mr and ns are, respectively, hr x ns and ns x ns mass 
matrices. Both right-handed neutrinos and singlet fermions have lepton number 
L = 1; consequently, the mass term S c /isS violates lepton number for two units. 

Introducing the left-handed basis N L = (u L , N£, S) t we have, after EW symme¬ 
try breaking, the following mass matrix 


M 


0 mjj 0 
m D 0 Mr , 
0 M^ ns ) 


( 2 . 2 ) 


with mo = vY v j\J 2. 

The mass matrix in eq. (2.2) can be diagonalized by means of the following 
unitary transformation 


U T MU = Ma => U\M ] M)U = M\ , N' L = WN l , (2.3) 


where M. a is the diagonal matrix referred to the mass eigenstates N' L . The first three 
eigenstates correspond to the standard light active neutrinos while the remaining 
nR + ns states are additional heavy sterile neutrinos. 

Following the standard seesaw calculation [52, 53] and assuming the hierarchy 
Mr » mp » Us, h is possible to extract the effective light neutrino mass matrix 

m u * m T D (Mj) * usM^rrio • (2.4) 

In the following we consider the case hr - ns = 3. Moreover, without loss of 
generality, we take Mr to be real and diagonal. We also work in a basis in which 
the mass matrix of charged SM leptons is diagonal. Within this framework, in the 
Us 0 limit the three light neutrinos are massless while the six heavy neutrinos 
can be recast into three pairs of Majorana particles with three (double degenerate) 
masses Mm, i = 1,2,3. On a general ground, from eq. (2.4) it follows that the 
order of magnitude of the light neutrino mass is m v ~ 0(/is x Assuming 

H,s ~ 0{l) keV, the model can accommodate sub-eV light neutrino masses with 
Y v ~ 0(1) couplings and Mr ~ (9(1-10) TeV seesaw scale. These order of magnitude 
estimates lie at the hearth of the ISS scenario. Small values of /is are expected by 
virtue of the’t Hooft naturalness criterium [54], since the limit ^s 0 increases the 
symmetry of the theory. Interestingly, the keV scale nicely fits the typical mass scale 
characterizing warm dark matter; in the ISS models with Ur + ns the spectrum - in 
addition to light and heavy neutrinos - also contains intermediate states with keV 
mass that are valuable warm dark matter candidates [55]. The estimates Y u ~ 0(1) 


4 




and M r ~ 0( 1 - 10) TeV represent the most relevant phenomenological properties 
of the ISS model since they allow for sizable (and, in principle, measurable) mixing 
effects between light active and heavy sterile neutrinos. This issue is particularly 
striking if compared with the typical high-scale characterizing the minimal type- 
1 seesaw [37-40] - that is Mr ~ 0( 10 15 ) GeV for order one Yukawa couplings - 
in which mixing effects, typically of order are negligible. As we shall 

explain in detail in section 3, the occurrence of both the peculiar conditions Y u ~ (9(1) 
and Mr ~ (9(1 - 10) TeV is of fundamental importance to determine the stability of 
the EW vacuum in the context of the ISS model. 

Of particular relevance for many phenomenological applications is the generalized 
Casas-Ibarra parametrization [56] 


Y v 


V2 

v 


V 'JMR vArt,:f/p M V, , 


(2.5) 


where \fm v is the diagonal matrix defined by the square roots of the eigenvalues 
corresponding to the three light neutrinos, 'rn vl with i = 1,2,3 hereafter, and VM is 
the diagonal matrix containing the square roots of the eigenvalues of M = Mrh^M^ 
whose diagonalization is defined by means of the transformation V T MV = M. R is 
an arbitrary 3x3 complex orthogonal matrix parametrized by three complex an¬ 
gles which encodes the remaining degrees of freedom. Finally, Dpmns corresponds 
to the unitary Pontecorvo-Maki-Nakagawa-Sakata (PMNS) leptonic mixing matrix. 
Assuming the standard picture with three neutrino flavors the matrix Dpmns can be 
parametrized as follows 


/ 


C12Y3 


S12C13 


Sise-iiop 


\ 


GpMNS - 


-S12C23 - C 12 S 13 S 2 3e iScp 
y S12S23 - Ci2Si 3 C23e i,5cp 


C12C23 - Si2Si3S2 3 e* ,5cp 
-S12S23 ” Sl2Sl3C2 3 e i<5cp 


C13S23 
C13C23 J 


( 2 . 6 ) 


where = cos % and s l3 = sin 9ij. In addition to the Dirac CP violation phase 
5 C p there are also two Majorana CP violation phases (not shown in eq. (2.6)). The 
latter are physical only if light neutrinos are Majorana particles, otherwise they can 
be always rotated away from the Lagrangian in the mass basis. In the following we 
omit the three phases of the PMNS matrix, since the Majorana phases are completely 
unknown and there are only preliminary hints about a non-zero Dirac phase. 


2.2 Bounds from low-energy neutrino data 

Global analysis of neutrino oscillation data, based on the latest results of the Daya 
Bay [57], RENO [58, 59], T2K [60, 61] and MINOS [62, 63] experiments, allowed 
to determine the oscillation parameters Am^, IA 777 . 3 -J (lAmjy, depending on the 
ordering), # 12 , 0 23 , 6* 13 with unprecedented high precision, thus opening the era of 
neutrino precision measurements. 
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In this work, we use the latest results of the Z/fit group [64], As customary, 
we define A m\- = m 2 ?: - m^, and we adopt the convention that results for the mass 
squared differences are reported with respect to the one with the largest absolute 
value. 


o Neutrino mass squared differences 

3-c C.L. ranges on the mass squared differences 


Au4/l(r 5 eV 2 = (7.02 -* 8.09) , 


f AmyiO-W 2 = (2.317 -* 2.607) NO 

I Am 2 2 /10- 3 * eV 2 - (-2.590 -2.307) IO 

(2.7) 

where the first (second) possibility refers to the assumption of normal (inverted) 
ordering. 


o Leptonic mixing matrix 

3-a C.L. ranges on the magnitude of the elements of the leptonic mixing matrix 
in eq. (2.6) 


sin 2 9 12 = (0.270 -> 

► 0.344) , 

(2.8) 

J (0.382 - 0.643) 

Sin 6 23 = i , Sin 013 = 

| (0.0186 ^ 0.0250) 

NO / x 
2.9 

[ (0.389 -* 0.644) 

\ (0.0188 -*■ 0.0251) 

IO 


o Unitarity 

In the basis where the charged lepton mass matrix is diagonal, the leptonic 
mixing matrix in the ISS model is given by the rectangular 3x9 sub-matrix 
corresponding to the first three rows of the matrix U defined in eq. (2.3), with 
the 3x3 block corresponding to the (non-unitary) ZTpmns- Bounds on the non- 
unitarity of the matrix ZTpmns were derived in [65-67] using an effective field 
theory approach. These bounds can be recast as follows 2 


^a/3 — 


E u°iU> 


2=4 


Pi 


|CpMNsCp 


PMNSl 


\$af3 (CpMNsCpjyjp^g)a/3| > (2.10) 

f (0.9979 -* 0.9998) < IO" 5 < 0.0021 ' 

< IO' 5 (0.9996 -> 1.0) < 0.0008 

< 0.0021 < 0.0008 (0.9947 ->■ 1.0) 


2 Strictly speaking, the bounds in eq. (2.10) are valid only if the masses of the sterile neutrinos 

lie above the EW scale (where they can be safely integrated out). This is always the case in our 

numerical analysis (see section 2.3). 
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o 


Additional constraints 


The absolute values of neutrino masses m ui are unknown. Cosmology sets the 
most stringent upper bounds using data from the Cosmic Microwave Back¬ 
ground (CMB) radiation, supernovae and galaxy clustering. Assuming the 
validity of the ACDM model [68], the Planck collaboration placed the upper 
bound Y,i m vi < 0.66 eV at 95% C.L. [68]; this bound becomes even more strin¬ 
gent adding data on the Baryon Acoustic Oscillation, Y,i m vi < 0.23 eV at 95% 
C.L. [68]. In our analysis we scan over the interval 10~ 4 eV m u i ^ 10 _1 eV for 
the mass of the lightest neutrino. 


2.3 Relevant parameter space and setup for the numerical analysis 
2.3.1 Target observables 

The presence of sterile neutrino states affects the SM charge current interaction via 
the mixing matrix in eq. (2.3). Going from gauge to mass eigenstates we have 


£ C c = — Yj vil + h.c. 

V ^ l=e,fi,r 


9 2 

V2 


9 


E Y^'fw-uu^ + h.c., 


( 2 . 11 ) 


where g 2 is the weak coupling constant. By means of these interactions, and de¬ 
pending on their masses and mixings with light active neutrinos, the presence of 
new sterile states can sizably affect numerous observables, like for instance leptonic 
and semi-leptonic decays (with a special focus on flavor violating processes) [69-78] , 
invisible Z boson decay width [79], Higgs boson decays [80-83], direct production in 
meson decay [84], 

In this paper, in order to investigate the stability of the EW vacuum in a region 
of the parameter space of particular interest for present and future experimental 
prospects, we focus on the lepton flavor violating process fi -> ey and the neutrino¬ 
less double beta decay (0u2f3 hereafter). 

As far as the radiative /i ey decay is concerned, the rate induced by the 
presence of sterile neutrinos is given by [85, 86] 


Br(|U -* ey) 


3ct 

327T 




i= 1 


m U 

Mw, 


( 2 . 12 ) 


where a is the electromagnetic fine structure constant and Mw the W mass. The loop 
function is given by G(x) = (10 -43x + 78x 2 -49x 3 + 4a; 4 + 18a; 3 lna;)/[3(-l + a;) 4 ]. The 
present experimental upper bound, reported by the MEG collaboration, is -*■ 
erf) < 5.7 x 10" 13 at 90% C.L. [87], 

The amplitude of the 0z/2/3 process is proportional to the so-called effective neu¬ 
trino mass, m v e Q. Current experiments (among others, GERDA [88], EXO-200 [89, 
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90], and KamLAND-ZEN [91]), put an upper limit in the range \m^\ < 140 - 700 
MeV. In the presence of sterile state the effective neutrino mass is given by [92] 


m 


eff 




m h 


i=1 


p 2 - ml- 


Y1 Uei m u i + p 2 [ f/g 4 


i— 1 


m v4 2 m j/5 

9 9 ^ eb 9 2 

P ~ m l 4 P ~ m u5 


+ U, 


m u 6 


e 6 o 9 

p 2 - ml 6 


+ f/. 


m u7 

el 9 9 

- m^ 7 


+ t/, 


"As 


eS p 2 -ml s 


+ U' t 


m v9 


e9 P 2 - mlg 


(2.13) 


where p 2 ^ -(125 MeV ) 2 is the momentum of the virtual neutrino. Notice that, since 
we are considering the regime m 2 i=4 9 »p 2 , heavy neutrinos decouple in eq. (2.13) 

and the dominant contribution to m u e ^ comes from the light active neutrinos. 


2.3.2 Strategy and first numerical results 

We perform a scan over the parameter space of the model, and our procedure goes 
as follows. First, we randomly generate i) the light neutrino masses m V i= 1 , 2,3 and 
the leptonic mixing angles $ 12 , # 23 , 6 13 according to the corresponding l-a intervals 
allowed by the analysis of present experimental data, ii) the entries of the matrices 
M r and ps in the intervals 10 2 GeV ^ M Ri ^ 10 2 TeV, 10 _1 keV ^ (ps)ij ^ 10 2 keV 
and in) the complex angles defining the arbitrary matrix R in the interval [0,27 t], 
Second, we reconstruct the full Yukawa matrix Y v using the generalized Casas-Ibarra 
parametrization in eq. (2.5). Finally, plugging back Mr , ps and Y v into eq. (2.2), 
we diagonalize the mass matrix J\4 in order to find the full 9x9 mixing matrix U. 
The phases of the mixing matrix are fixed using eq. (2.3), by means of the condition 
m„i ^ 0 for all i. As a consistency test, for each point of the scan we check that 
the mass matrix A4, randomly generated as discussed above, correctly reproduces 
after diagonalization light neutrino masses and mixing angles in agreement with the 
bounds discussed in section 2.2. Equipped by these results we can easily compute the 
branching ratio Br (p -*■ eq) in eq. (2.12) and the effective neutrino mass in eq. (2.13). 
We show our results in fig. 1 for the normal ordering and in fig. 2 for the inverted 
ordering. In the upper (lower) panels of both figures we show the branching ratio 
Br(/i - 2 - eq) (the effective neutrino mass m[(|) as a function of the lightest heavy 
neutrino mass m u4 (plot on the left) and the trace of the Yukawa couplings Tr(YjY u ) 
(plot on the right). Few comments are in order. 

1. Normal and inverted ordering produce very similar distributions considering 
the radiative decay Br (p eq). This is caused by the well-known fact that the 
contribution of light active neutrinos is strongly suppressed by the extremely 
small value of light neutrino masses. In the ISS model a non-zero contribution 
to Br(p eq) is entirely generated by the additional heavy neutrinos, and 
controlled by the mixings 11^,11^ ~ mo/M r (see eq. (3.11) in section 3.2). We 
notice that in our scan we can obtain a signal close to the present experimental 
bound even considering m v4 as large as 10 TeV and Tr(YjV,) as small as 
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> 
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Figure 1: Upper panel. Branching ratio for the decay process n ey as a function of 
the mass of the lightest sterile neutrino (left panel) and the trace of the Yukawa cou¬ 
plings Ti :(YjY v ) (right panel). The blue horizontal line represents the upper bound set 
by the MEG collaboration [87]. Lower panel. Effective neutrino mass as a function 
of the mass of the lightest sterile neutrino (left panel) and the trace of the Yukawa 
couplings T t(YJY v ) (right panel). The blue solid (dashed) line represents the upper 
bound (future sensitivity) of the EXO-200 experiment [90], Details about the nu¬ 
merical scan are given in the text. All points comply with the bounds discussed in 
section 2.2. 


10 3 . For completeness we show in the left panel of fig. 3 the result of our 
numerical scan in the plane [Tr (Y^Y u ),m v f\. We mark in dark cyan points with 
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- 
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- 


10 3 10 4 10 5 10~ 5 10" 4 10“ 3 10“ 2 KT 1 10 ° 10 1 10 2 

m v4 [GeV] TRFjTy) 


Figure 2: Same as in fig. 1, but considering the inverted ordering. 


Br (n -> 07 ) ^ 10 -13 . Points with large Yukawa couplings (e.g. Tr (YjY' u ) > 0.5) 
and sizable /a -*■ 07 rate are generated in the whole interval of analyzed masses 
for the right-handed neutrinos. 

2. Normal and inverted ordering produce completely different distributions con¬ 
sidering the effective neutrino mass mfig. In this case the contributions of 
additional heavy neutrinos decouple since their masses are much larger than 
the typical momentum scale p 2 ^ -(125 MeV) 2 . Therefore, in our numerical 
scan the ISS model resembles the typical scenario with only three light active 
neutrinos. The situation is well represented by the right panel of fig. 3 where 
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we show the effective neutrino mass for the normal and inverted ordering as a 
function of minimal neutrino mass. The normal ordering is suppressed since 
the largest neutrino mass is multiplied by the small value of S 13 . However, in 
both cases the effective neutrino mass is close to the future sensitivity of the 
EXO-200 experiment. 



Tr (Y}Y V ) 


m v 1 [eV] 


Figure 3: Left panel. Distribution of our numerical scan in the plane 
[Tr(YjY v ), m^] considering the normal ordering (the inverted ordering gives an anal¬ 
ogous result). All points comply with the bounds discussed in section 2.2. Points with 
Bi :(/i -»• ey) ^ 10 ~ 13 are marked in dark cyan. On the qualitative level, the plot shows 
that Yukawas Tr (YjY v ) > 0.5 arise in the whole range of analyzed masses for the 
right-handed neutrinos. Right panel. Effective neutrino mass as a function of the 
lightest neutrino mass considering both normal and inverted ordering. Our scan cor¬ 
rectly reproduces, as expected if the heavy neutrinos decouple if compared with the 
typical virtual momentum p 2 - -(125 MeV) 2 , the well-known result characterizing 
the presence of only three light active neutrinos. 


2.3.3 Additional remarks 

Let us close this section summarizing further predictions and constraints on the in¬ 
verse seesaw scenario. The aim of the following discussion is to provide additional 
motivations enforcing the phenomenological relevance of EW-scale right-handed neu¬ 
trinos with 0(1) Yukawa couplings. 

Collider searches at the LHC. At the LHC right-handed neutrinos with a mass not 
far above or below the Higgs mass and with a sizable Yukawa coupling Y u > 0( 10“ 2 ) 
affect the Higgs decay h -* llvv (see [80] for a recent analysis). Present bounds 
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hold in the range 60 ^ Mr ^ 200 GeV with 10~ 2 2. 3 Furthermore, the CMS 

collaboration placed upper limits on the active-sterile neutrino mixings in the same 
mass range for M R considering direct production of heavy neutrinos [93, 94]. 4 As 
stated in the introduction, these values of masses and couplings may have an impact 
on the stability of the EW vacuum, thus providing an additional motivation for the 
analysis that we shall perform in the next section. 

Fit of LEP data via oblique parameters. The fit of LEP data still provides today 
an important constraint on beyond the SM physics. The presence of additional sterile 
neutrinos modifies the oblique radiative corrections [96-103]. In [79] it was shown 
that right-handed neutrino masses of the order M Ri ~ (9(10) TeV, together with 
violation of unitarity of the order e a p ~ 10 -3 -f-10~ 6 , can improve the fit of LEP data 
with respect to the SM. 

Leptogenesis. In the inverse seesaw scenario the decay of (nearly degenerate) 
heavy Majorana neutrinos can realize the so-called resonant leptogenesis [104]. Re¬ 
markably, resonant leptogenesis can be realized with heavy Majorana neutrinos as 
light as 1 TeV [117] (in contrast with the usual thermal leptogenesis, realized in the 
type-I seesaw, in which M n > 10 9 GeV [106]). 

Naturalness. On a general ground, whenever a threshold with particles of mass 
M coupled to the Higgs with strength £ is present, quantum corrections generate a 
contribution 5m 2 H » £ 2 M 2 /167t 2 to the renormalized Higgs boson mass. If 5m 2 H » v 2 , 
an unnatural cancellation between 8m 2 H and the bare Higgs mass is required in order 
to reproduce the observed value of the Higgs boson mass. The condition 5m 2 H < v 2 
can be used as a criterium to construct natural model of new physics [107, 108]. In 
the context of inverse seesaw models, the scale M is the mass of right-handed heavy 
neutrinos Mr, and the coupling £ is the Yukawa coupling Y v . TeV-scale values of 
M r and (9(1) Yukawa couplings satisfy the naturalness condition. 

To sum up, the presence of heavy neutrinos with a mass not far above the 
EW scale and a sizable Yukawa coupling with the Higgs boson has extremely rich 
phenomenological consequences. Motivated by the results achieved in this section, 
we are now in the position to tackle the second part of our analysis in which we aim 
to investigate the impact of heavy neutrinos on the stability of the EW vacuum. 

3 Stability of the electroweak vacuum and low-scale seesaw 

The key ingredient in the study of the stability of the EW potential is represented by 
the quantum effective potential V e Q. In particular, we search for the instability scale 
A corresponding to the Higgs field value for which the potential becomes smaller 

3 These bound were obtained in [80] considering a simplified setup with only one light flavor of 
heavy neutrinos. Consequently, here Y v indicates the corresponding Yukawa coupling to the Higgs 
doublet. 

4 See [95] for prospects at future lepton colliders. 
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than its value at the EW minimum. If such scale does not exist the vacuum is 
called absolute stable , otherwise the Higgs vacuum is not the global minimum, and 
a quantum tunnelling to the true vacuum may occur. The decay probability can 
be computed from the bounce solution of the euclidean equations of motion of the 
Higgs field [18, 19]. If the corresponding lifetime is much bigger that the age of the 
Universe the EW vacuum is called metastable. This can be translated into a lower 
bound on the Higgs effective quartic coupling A e s, driven to negative values, that at 
leading order reads as 




8tt 2 1 
3 log(r^i) ' 


(3.1) 


where r is the age of the Universe r = 4.35 x 10 17 sec [109] and /a is the renormalization 
scale of the RG running. The Higgs effective quartic coupling will be defined in the 
next section. The SM quantum corrections at zero temperature to the metastable 
condition have been computed in [110] but are not considered in this work. The 
measured values of the Higgs and the top masses place the SM in a metastable 
position in the phase space diagram [111]. The inclusion of extra fermionic degrees 
of freedom, as right-handed neutrinos in seesaw extensions, can clearly change this 
picture possibly driving the model to an instability region. Therefore, by requiring 
the metastability of the EW vacuum, together with the perturbativity of the gauge, 
scalar and Yukawa couplings up to the Planck mass, we can constrain the parameter 
space of low-scale seesaw models. 

In the following sections we will describe the theoretical tools needed in our sta¬ 
bility analysis: the Higgs effective quartic couplings which determines the instability 
scale, the RGEs in the MS describing the running of the couplings from the EW 
scale up to the Planck scale and, finally, the matching conditions which provide the 
initial values of the MS parameters from the EW physical observables. 


3.1 The Higgs effective quartic coupling 

From a preliminary analysis one can show that the vacuum instability appears at a 
scale much bigger than the EW minimum. This allows us to consider the effective 
potential at large field values, (j) » v, and to use the following approximation 

(3.2) 

where the Higgs quadratic mass term has been neglected and t is the logarithm of the 
renormalization running scale /a. The effective quartic coupling A e fr is extracted from 
the RG-improved effective potential at two-loop order in the SM [112] and at one- 
loop for the right-handed neutrino corrections, computed in the MS renormalization 
scheme and in the Landau gauge. The complete two-loop expression is too lengthy 
to be given here. We show, instead, the effective quartic coupling at one-loop order 
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in the SM [113] 


A« e 4r(<) < A(f) 


167T 2 


E N P K l(t) 


, K p (t)e 2r ^(j) 2 ^ 

1 “TM 1 p 


p=W,Z,h,x,t 

where the ^coefficients are summarized in table 1 and T(t) is defined as 

m= Cdt'^f) , 

Jo 

with 7 the anomalous dimension of the Higgs held. 



t 
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Z 

h 
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X° 
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1 
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3/2 

hip 
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(«£ + #)/■ 4 

3A 

A 

A 


(3.3) 


(3.4) 


Table 1: SM p-coefficients entering in eq. (3.3). 


The contribution of the heavy neutrinos to the RG-improved effective potential 
is [114] 


AW.i) = -—J £ 

6Z7T i=l,2,3 


M 2 - 

In 


(<M) 


p(ty 


(3.5) 


where 9 m (t) = 6{t -In M Ri /p 0 ) and M. vi (<f),i) are the three (double degenerate) non¬ 
zero eigenvalues of the 9x9 mass matrix 


M4t) 


0 W(f) T 0(f)/V2 0 ' 

Y v (t)^(t)/y /‘.2 0 , 

0 Mr?, 0 


(3.6) 


with Mr? = diag(MRi, Mr 2 , Mr 3 ) and 0(f) = e r b)0. Compared with eq. (2.2), we are 
considering the /xg -► 0 limit in which the three light neutrinos are massless. 

The factor of two in eq. (3.5) comes from the fact that each non-zero eigenvalue 
is double degenerate. 

The RG running parameter t can be chosen in such a way that the convergence 
of perturbation theory - otherwise spoiled by the presence of large logs - is improved. 
We follow the prescription /x(t) = 0 widely used in vacuum stability analyses. 

The contribution of heavy neutrinos to the effective potential produces two dis¬ 
tinctive effects on the RG-evolution of the effective quartic coupling in eq. (3.3). 
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Above each threshold M® and at high values of renormalization scale /i(t) » 
M Ri it contributes explicitly to A e ff introducing the correction 


A&(<M) 


167T 2 


4 >4i(t) 


In 


K^(t)e 2r h)0 2 

W ) 2 



(3.7) 


where K V i(t) is the coefficient of the ^-dependent part of M U i((p,t). 


o Below each threshold M Ri the corresponding heavy neutrinos are integrated 
out. The matching produces the threshold correction to the effective potential 




32t r 2 


2 MUM 


M 2 - 

In 


om) 


Mfu 


(3.8) 


which translates into a threshold correction AA{£ to the Higgs quartic coupling 
at the M Ri mass scale which can be extracted from the 0 4 term of eq. 3.8 and 
it is explicitly given by 


aa 


1 dAV£(<j>,t) 
6 d(p 


(3.9) 


3.2 The matching conditions 

The RG equations employed in this work are computed in the MS renormalization 
scheme and must be equipped with suitable initial conditions for the running pa¬ 
rameters evaluated in the same scheme. These parameters are expressed in terms 
of physical observables, defined in the on-shell OS scheme, through appropriate re¬ 
lations called matching conditions. For a generic parameter a the matching at the 
scale n is given by 


c*(aO = «os - <5«os + , (3.10) 

where a(/i) and «os are respectively the MS and the OS parameters, while tfayjg 
and 5«os are the corresponding counterterms in the two schemes. The difference 
between the two of them is an ultraviolet-finite correction. 

The SM parameters obtained from the matching procedure are (X,yt, ( 71 ,( 72 ), 
while the QCD coupling constant g 3 is directly extracted from the a^(Mz) which is 
already defined in the MS scheme. The physical observables used to compute the 
input values of the MS parameters are given in table 2. 

The details of the strategy can be found in [111] where the SM two-loop (NNLO) 
corrections to the matching conditions have been discussed. In particular, in [111] 
a complete two-loop analysis has been performed in the EW sector and the N3LO 
(three-loop) pure QCD effect has been included in the matching of the top Yukawa 
coupling and the strong coupling constant. The running of the latter from Mz to 
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Name 

Value 

Description 

M\y 

80.384 GeV 

W boson pole mass 

M z 

91.1876 GeV 

Z boson pole mass 

M h 

125.09 GeV 

Higgs boson pole mass 

M t 

173.34 GeV 

Top quark pole mass 

v e (v/^r 1 / 2 

246.21971 GeV 

Higgs vev from the y decay 

az(Mz) 

0.1184 

MS QCD structure constant (5 flavors) 


Table 2: Physical observables used to extract the SM parameters in the MS scheme 
through the matching procedure. For the Higgs mass we used the latest result [3], for 
all the other parameters we refer to [111]. 


M t , which is the starting scale of our stability analysis, has been performed including 
the QCD four-loop /3 function. 

In low-scale seesaw extensions, the right-handed neutrino corrections to the 
matching conditions can be important and must be taken into account. In our 
analysis we have considered all the SM results given in [111], supplemented by new 
physics contributions computed at one-loop order from the Lagrangian in eq. (2.1). 
In particular, the additional neutrinos introduce corrections to the masses of the 
gauge bosons M z , M w , the Higgs M h , the quark top M t and to the muon decay, 
from which the Fermi constant G^ is extracted. These corrections depend on the 
masses of the heavy neutrinos, their interactions with the SM fields mediated by the 
Yukawa couplings Y u (obtained from the parameterization in eq. (2.5)), and the full 
9x9 mixing matrix U. Due to the mass hierarchy Mr » mo, and to the smallness 
of the light neutrino masses - which can be safely neglected in the computation of 
the matching conditions at the EW scale - the mixing matrix U can be expanded in 
the ratio mo/M r, with ps 0, as 


U 


1 
0 

-Mjfrrio 


,/2 n ^ D Mji 

i 

n/2 
1 

n/2 


73 < M n 

_ i_ 

n/2 

i 

V2 


+ 01 


m 


D 


,K, 


(3.11) 


Notice that the PMNS block has been set to the unit matrix, consistently with the 
approximation ^ 0. 

Moreover, we have verified the Appelquist-Carazzone theorem in the new physics 
sector. In particular, we have checked the decoupling of right-handed neutrino con¬ 
tributions from the matching conditions in the limit of large Majorana masses. This 
has to be expected since these masses are not generated by spontaneous symmetry 
breaking of the Higgs field. 
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Concerning the matching conditions of the Yukawa matrix Y u , computing pertur¬ 
bative corrections to the matching conditions at the EW scale will not considerably 
improve the precision on the determination of Y v in the MS scheme (since in any case 
Y v turns out to be related to unknown parameters by means of eq. (2.5)). Therefore, 
for the sake of simplicity, we decided to match the MS Yukawa matrix Y v to its OS 
version at the tree-level, namely Y u (M t ) ^ Y v . 


3.3 The renormalization group equations 

All the dimensionless couplings (X,y t ,gi,Y u ) are evolved from the top-mass scale, 
M t , up to the Planck scale using the three-loop (NNLO) RG equations for the SM 
parameters and the two-loop (NLO) /3 functions for the Yukawa matrix Y v . Here 
Qi stands for the three gauge coupling constants and we have retained only the top- 
quark contribution in the SM Yukawa sector. The system of coupled RG equations 
is then given by 


d\(t) 

dt 

dyt(t) 

dt 

dgijt) 
dt 

dY v {t) 

dt 


P\(\yt,g%,Yu) , 

(3.12) 

Pyt (5 Ut 1 gi: 4 V ) 5 

(3.13) 

Pgii^iVti giiYj) i 

(3.14) 

pY v ( x xyt,gi,Y») , 

(3.15) 


where the /3 functions are computed in perturbation theory in the MS renormalization 
scheme. Due to their lengthy expressions, we present only the one-loop corrections 
to the r.h.s. of eqs. (3.12-3.15), namely 


Px = « [24A 2 + A (l2 y 2 + 4Tr (YjY v ) - - 9<? 2 2 ) - 6 yf - 2Tr (YjYj' 


27 4 9 4 
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+ 200 ^ + 8 92 + 20 9i92 


Pvt = K 


~3i - \al - 8sl + \y 2 t + TrK'FJ 
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Pgi = K—gi , Pg 2 = -K—gl , Pg 3 = -Klgl , 


Py„ = K 


r/__9_ 2 _ 9 
A 20^ 4 


gl + 3 y 2 + Tr 



V ) 


v;. i -i;.!;'!':. 


(3.16) 

(3.17) 

(3.18) 

(3.19) 


where k = 1 / (167T 2 ) and the abelian gauge coupling is given in the GUT normal¬ 
ization gi = y5/3 gy- Notice that the RG equations given above are defined for a 
renormalization scale fi bigger than any particle mass of the model. For lower scales, 
heavy degrees of freedom are integrated out and the corresponding coupling must be 
removed by hand from the /3 functions. Indeed, in the MS renormalization scheme 
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the decoupling of heavy degrees of freedom is not automatic and has to be explicitly 
implemented at the different particle thresholds. 

As far as the vacuum stability is concerned, the right-handed neutrinos behave 
like the top-quark and drive A to negative values faster than the SM case. This is 
clear from eq. (3.16). On the other hand, their impact on the top Yukawa coupling - 
see eq. (3.17) - is to increase y t all along the RG evolution. Actually, this is another 
source of vacuum destabilization with respect to the SM picture, because a bigger 
value of y t has a bigger decreasing effect on the Higgs quartic coupling. Nevertheless, 
the overall behavior of the top-quark Yukawa coupling is dictated by the large and 
negative QCD corrections which lead to a decreasing y t . These features can be easily 
deduced from fig. 4, left panel, where the running of the SM couplings is depicted. 
Here solid lines represent the evolution in the seesaw extended scenario, while dashed 
curves correspond to the pure SM. 

Contrary to the top Yukawa case, the QCD contributions are obviously absent - 
at least in the leading one-loop approximation - from the evolution of Y u , and the 
Tt(yJy u ) term, which affects /3\ e (3 yt , is always increasing. This feature, shown in 
the right panel of fig. 4, has a negative impact both on the vacuum stability and 
on the perturbativity of Y u which can be violated, if \Y vi j\ > \/47r, during the RG 
evolution for sufficiently big values of the Yukawa coupling at the EW scale. In the 
same figure the decoupling of the heavy right-handed neutrinos below their mass 
thresholds is also manifest. Indeed, for /i « M Til , the N l R neutrino is integrated out 
and does not contribute to the RG running: the corresponding row in the Yukawa 
matrix Y v is frozen, and enters in the [3 functions only above the threshold scale 
Mm as shown in eq. (3.20) where we mark with a generic x non-zero entries for the 
Yukawa matrix Y u . 

SM EFTi EFT 2 Full theory 
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Finally, we show in fig. 5 the evolution of the effective quartic coupling A e ff in two 
different scenarios. In the left panel the inverse seesaw is realized with Tr (YJY V ) - 
0.36, while in the right panel the Yukawa matrix is such that Ti :(YjY v ) - 0.6. In the 
latter case the effects of Y v , which affects the RG running above the threshold scales, 
are quite large and A e fr is driven outside the metastability region below the Planck 
scale. 
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Figure 4: Left panel. RG evolution of the SM couplings. Solid lines take into account 
the effects of right-handed neutrinos in the inverse seesaw model with Tr (YjY v ) - 
0.36 at the EW scale. Dashed lines represent the running of couplings in the SM. 
Right panel. Evolution of Tr(YjY u ). The heavy right-handed neutrino thresholds are 
explicitly shown. 



RGE scale y [GeV] 



RGE scale // [GeV] 


Figure 5: RG evolution of the effective Higgs quartic coupling in the inverse seesaw 
model in two different setup: T t(yJy v ) - 0.36 (left panel) and Ti :(YjY u ) - 0.6 (right 
panel). In the latter case the Yukawa couplings have a sizable impact on A e g- and the 
metastability bound is violated below the Planck scale. 
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4 Results and discussion 


In this section we discuss the results of our analysis. First, let us briefly remind 
our strategy. The Yukawa matrices generated following the prescription outlined 
in section 2.3 enter as initial conditions, together with all the other SM external 
parameters, in the solution of the RGEs in eqs. (3.12-3.15). From the effective 
potential, improved by the running couplings previously computed, we extract the 
Higgs effective quartic coupling A e fr in eq. (3.3). Finally, we use eq. (3.1) to check 
whether the analyzed points violate the metastability bound on the EW vacuum. 

We show our results in fig. 6 for the normal ordering and in fig. 7 for the inverted 
ordering. In order to make contact with phenomenology, in both cases we present 
the impact of the metastability bound with respect to the observables targeted in 
section 2.3.1, namely the branching ratio Br(/i -»■ eq) (left panel) and the effective 
neutrino mass (right panel). The red points are excluded by the metastability bound: 
the EW vacuum would decay too fast in the true vacuum of the EW potential. 


Inverse seesaw: Normal Ordering 



Tr(F v W v ) Tr(FjF v ) 


Figure 6: Impact of the vacuum stability analysis on the branching ratio Br(/i -> eq) 
(left panel) and the effective neutrino mass (right panel) in the ISS model. All points 
comply with the bounds discussed in section 2.2. Red points violate the metastability 
bound discussed in section 3. 

Our numerical analysis clearly indicates that points with Yukawa couplings such 
that Tr (YjY u ) > 0.4 are excluded. This bound does not depend on the assumed 
hierarchy, since light neutrino masses are irrelevant. Most importantly, the excluded 
points lie in a region of the parameter space that is close to the present bounds and 
future experimental sensitivities for both the analyzed observables. Moreover, as 
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clear from the left panel of fig. 3, the metastability bound applies in the whole range 
of analyzed masses for the right-handed neutrinos. This is an interesting piece of 
information since, for instance, searches for heavy neutrinos with m .,,4 ~ (9(100) GeV 
and sizable Yukawa couplings are currently ongoing at the LHC (see section 2.3.3). 


Inverse seesaw: Inverted Ordering 



ltr 5 icr 4 icr 3 to - 2 icr 1 io° to 1 io 2 icr 5 10- 4 icr 3 icr 2 io-‘ io° to 1 io 2 

Tr(K v + r v ) Tr(Y f v Y v ) 

Figure 7: Same as in fig. 6, but considering the inverted ordering. 


Before moving to conclusions, we refer the reader to the appendix A for a gen¬ 
eralization of these results to the linear and double seesaw models. In fact, one 
might be concerned about the model-dependence of our bound, since we worked in a 
well-defined model. By extending our analysis to the two aforementioned models, we 
intend to prove that the significance of the metastability bound is a general aspect 
of low-scale seesaw models. 

5 Conclusions and outlook 

In this paper we analyzed the stability of the EW vacuum in the presence of low-scale 
seesaw models. For definiteness, in the main part of this manuscript we focused on 
the ISS model, but in appendix A we generalized our results to other popular low- 
scale scenarios. In the ISS model lepton number is broken at the scale ns, and the 
scale of light neutrino masses is given by m u « fis( v Y u /M r ) 2 where Mr is the mass 
of the heavy neutrinos. Our analysis can be divided in two parts. 

First, in section 2, we carefully chose the parameter space of our numerical 
scan in order to comply with low-energy neutrino data. We performed this selec¬ 
tion by means of the Casas-Ibarra parametrization, and we considered the intervals 
10 2 GeV ^ Mr ^ 10 2 TeV and 10” 1 keV ^ ns ^ 10 2 keV in order to generate order one 
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Yukawa couplings. As expected, the output of our numerical scan shows that the 
ISS model presents an extremely rich low-energy phenomenology, and we decided to 
focus our attention on the lepton flavor radiative decay fi -»• ey and on the 0n2 (3. 

Second, in section 3, we computed the contribution of the heavy neutrinos to the 
RG evolution of the dimensionless SM couplings and to the RG-improved effective 
potential. On the technical level, we worked including the most updated results for 
the computation of (3 functions and matching conditions in the SM, and we added 
the contributions of the right-handed neutrinos at two-loops in the (3 functions and 
one-loop in the matching. 

With such heavy artillery, we analyzed the Yukawa couplings generated in sec¬ 
tion 2 by imposing the metastability condition on the lifet im e of the EW vacuum. 
Our conclusions are summarized in fig. 6 and fig. 7: Yukawa couplings such that 
Tr (YjY u ) > 0.4 are excluded by the metastability condition. The bound applies 
in the whole range of analyzed right neutrino masses, and affects a region of the 
parameter space close to present or future sensitivities for both -> ey rate and 
0v2(3. 

To sum up, we argue that the metastability bound represents an important 
consistency condition that should be included in all the phenomenological analysis 
of low-scale seesaw models featuring EW-scale right-handed neutrinos with (9(1) 
Yukawa couplings. 
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A Linear and double seesaw 


In full generality, we introduce a Majorana mass term for the right-handed neutrinos 
and a Yukawa coupling for the fermionic singlets 


£ 


iN R ry li (d li N R )+iS'y li (d li S) 

NrY u H^L + SCYsH^L + N^M r S + l -NgM N N R + + h.c. 


• (A.l) 
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In the left-handed basis Nl = 
the following mass matrix 


(yr, Nr, S) t we have, after EW symmetry breaking, 


M 


0 m/) rrig 

uid Mn Mr 


K m s Mr p 5 t 


(A.2) 


with m-D = vY u /\/2 and ms = vYs/V2. 

o The linear seesaw [115] corresponds to /is = 0, M N = 0 


M 


LSS 


0 m^ ^ 

m D 0 M r 

ms Mr 0 


(A.3) 


In this case the only source of lepton number violation comes from ?%. Fol¬ 
lowing the standard seesaw approximation we find 

m u « + m^M^m D . (A.4) 

Generalizing the Casas-Ibarra parametrization [56], we obtain [53] 

V 2 

Y v = —M R (mf) -1 t/p MNS Vr^Ay^t/^ MNS , (A.5) 

where A is a 3 x 3 matrix satisfying the equation A + A T = 1. Consequently, 
An — 1/2 and Aij — —Aj 


o The double seesaw corresponds to ms = 0 

/ 


M 


DSS - 


0 777./] 


mr, Ain Mr 

0 MJ /< : ) 


(A.6) 


The Majorana mass term Mat, in addition to /is, violates lepton number for 
two units. We assume the hierarchy Mn » Mr » mr> » /is- Integrating out 
the heavy fields Nr, we find the effective Lagrangian 

■Cdss = (A.7) 

+ l -J^H'(Y^M N 'M R ) v S, + isf - I§<%S + fc.c. . 

In the left-handed basis AA = ( u L ,S) T eq. (A.7) corresponds, after EW sym¬ 
metry breaking, to the mass matrix 


l-m T D M^m D -m T D M^M R \ 
\ -MlMfmo Ps - MIM-JMr ) 


23 








After block-cliagonalization we find 


V T M$ SS V 


m v 0 

0 Af^gayy 


where 


(A.9) 


m u = -m^M N 1 m D -(m, r ^M N 1 M R )(^ s -M R M N 1 M R ) 1 {M R M^m D ) , (A.10) 

A/heavy = f^S ~ M R M N Mr ■ (A. 11) 

From eq. (A. 10) it follows that m v -> 0 if /is -»• 0 since the type-1 contribu¬ 
tion mJjM^}rrio cancels out between the two remaining terms. Considering a 
perturbative expansion in fj,g, we find 

m v « m'i ) (M R )~ 1 fj lS M R -m D . (A.12) 


As a result, the Casas-Ibarra parametrization is analogous to eq. (2.5). The 
unitary mixing matrix V in eq. (A.9) has the general structure 


V 


Vl^BB t B 

-Rt Vl - B^B 


(A.13) 


and, at the lowest order, we find B* = m^(AfJ) -1 . For the sake of simplicity, we 
focus on the case with n R = n$ = 3, and we assume M R = diag(Af ffi ), and M N = 
diag(Afjvi) with i = 1,2,3. In this case the heavy block in eq. (A.9) simplifies to 
AfAf;^ 1 M R = diag(Af^ i /Afjvj). In eq. (A.7) the interactions between the lepton 
doublets Lj and the three singlet fermions S% are mediated by an effective 
Yukawa matrix Y v = Mj, M^} Y v . The previous assumption implies (Yv)ij = 
Mm ( Y v ) jj / M Nl . 


A.l Linear seesaw: results 

For definiteness, we focus on the case with n R - n R - 3. Considering the low-energy 
neutrino data, we follow the same strategy outlined in section 2 (see in particular 
section 2.3). We make use of the Casas-Ibarra parametrization in eq. (A.5) to numer¬ 
ically reconstruct the Yukawa matrix Y v , and we randomly scan over the intervals 
10 2 GeV ^ ^ 10 2 TeV, 10 _2 keV {jns)ij ^ 10 2 keV, and 10 _1 ^ A l3 ^ 10 2 . As 

done for the inverse seesaw, we discard points unable to comply with the bounds 
discussed in section 2.2. As far as the stability of the EW vacuum is concerned, in 
the ms -*• 0 limit the mass matrix A4 lss hi eq. (A.3) reduces to the same structure 
already studied in the inverse seesaw case (see eq. (2.2)). Consequently, in the def¬ 
inition of the Higgs effective quartic coupling we employ the same RG equations, 
matching conditions and effective potential used in the inverse seesaw analysis. We 


24 




Linear seesaw: Normal Ordering 
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Tr(FjF v ) = 0.4 


Figure 8: The same as in fig. 6 but considering the linear seesaw. 

show the final result of our analysis in fig. 8 considering a normal ordering of light 
neutrino masses. 

As expected, we find the same quantitative conclusion if compared with the 
inverse seesaw case. Yukawa couplings such that Tr(YjY„) > 0.4 are excluded by 
the metastability bound. Remarkably, considering both the lepton flavor violating 
process /i -»■ eq and the 0z^2/3, these points lie in a region of the parameter space 
close to present or future sensitivities. Therefore, we conclude that the metastability 
bound represents an important consistency condition that should be included in all 
the phenomenological analysis of the linear seesaw model featuring EW-scale right- 
handed neutrinos with 0(1) Yukawa couplings. 

A. 2 Double seesaw: results 

We follow the same approach already exploited for the inverse and linear seesaw 
models. However, in the double seesaw case there are few remarkable differences. 
We use the Casas-lbarra parametrization in eq. (2.5) to sample the Yukawa matrix 
Y u , and we randomly scan over the intervals lOTeV ^ ^ 10 3 TeV, 10 9 GeV ^ 

M m ^ 10 11 GeV, and lOkeV ^ (^s)ij ^ 10 4 keV. Few comments are in order. First, 
notice that this choice of parameters - optimized in order to obtain 0(1) Yukawa 
couplings - respects the hierarchy M N » M R » m D » fig assumed above (see 
discussion below eq. (A.6)). Second, we expect the following order of magnitude 
estimates: for the mixing parameter in eq. (A.13), B ~ (9(10~ 3 ); for the mass of 
the heavy neutrinos in eq. (A. 11), Mh eavy ~ 0(1) GeV; for the effective Yukawa 
coupling in eq. (A.7), Y v ~ 10 _5 Yj,. Armed with these numbers, we can outline as 
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follows. At large renormalization scale values, /i » M Nll the model is described 
by the full Lagrangian in eq. (A.l) (with Yg = 0). In terms of the RG running, 
the only relevant parameter is the Yukawa matrix Y u describing the interactions 
between the Higgs doublet, the lepton doublets and the three right-handed heavy 
neutrinos. At this stage, the situation is formally equivalent to the familiar type-I 
seesaw. 0 There is, however, one remarkable difference. In the type-I seesaw right- 
handed neutrino masses M Ni ~ (9(10 9 ) GeV require, in order to reproduce low- 
energy neutrino phenomenology, small Yukawa couplings (typically Y u ~ 10 -5 ). As a 
consequence, the impact of the interactions NgY u H^L + h.c. on the running of the 
Higgs quartic couplings is negligible. In the double seesaw case, on the contrary, 
we are allowed to consider (9(1) Yukawa couplings since the mass of light neutrinos 
is set by fig (see eq. (A. 12)). Below the thresholds M Ni the heavy right-handed 
neutrinos are integrated out, and eventually the model is described by the effective 
Lagrangian in eq. (A. 7). Given our choice of parameters, in this region the running is 
approximately equivalent to the pure SM since for the effective Yukawa interactions 
~LPH*Y V S we expect Y v ~ 10" 5 Y v . 

We summarize our results in fig. 9. In the left panel, we show the running of the 
effective Higgs quartic coupling for a specific realization of the double seesaw model 
with Tr (YjY u ) - 0.62. Above the thresholds the Yukawa couplings Y v sizably 
affect the running of A e fr eventually violating the metastability bound before the 
Planck scale. In the right panel we show the result of our numerical scan focusing on 
the effective neutrino mass relevant for the 0z/2/3. The most striking difference with 
respect to the inverse and linear seesaw models (see, respectively, figs. 6, 8) is that 
the presence of additional neutrinos with mass Mh ea vy ~ 0(1) GeV gives a sizable 
contribution to As a result, numerous points in our numerical analysis are close 
to (or even exceed) the present experimental bound. We find that Yukawa couplings 
such that Tr (YjY u ) > 0.6 are excluded by the metastability bound. 

Before concluding, let us stress that the aim of this discussion is not to be 
exhaustive since a careful phenomenological analysis is clearly far beyond the scope 
of this appendix. The most important message of our discussion is that even in the 
presence of very heavy right-handed neutrinos (with mass of the same size as the one 
expected in standard type-I scenarii) large modifications of the effective Higgs quartic 
couplings are a concrete and realistic possibility. As a final remark, notice that double 
seesaw models represent an appealing setup to realize leptogenesis [116, 117]. With 
this respect, a more detailed study of the parameter space will be addressed in a 
forthcoming publication. 

5 In concrete, the contribution of each heavy right-handed neutrino to the effective quartic cou¬ 
pling is given by eq. (3.7) (divided by two, since now there is no double degeneracy) while for the 
/3 functions we exploit the same two-loop expression already discussed in section 3.3 (see eq. (3.19) 
for the one-loop approximation). 
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RGE scale p [GeV] 


Tr (YX) 


Figure 9: Left panel. Running of the effective Higgs quartic coupling in the double 
seesaw model. Right-handed neutrinos with mass Mm ~ 0(1O 9 ) GeV (vertical dashed 
lines) sizably affect the RG evolution thanks to 0{l) Yukawa couplings. In the an¬ 
alyzed case, Yukawa couplings such that Tr(YjV„) 0.62 exceed the metastability 
bound below the Planck scale. Right panel. Effective neutrino mass as a function of 
the trace of the Yukawa couplings, Tr (YjY„). 
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